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Purpose: The objective of this study is to evaluate the effects of heat and massage 
application on autonomic nervous system. Materials and Methods: One hundred 
thirty-nine subjects volunteered and completed this study. Heat and massage was 
daily applied for 40 minutes, 5 days a week for 2 weeks. Primary-dependent mea- 
sures included heart rate variability, sympathetic skin response, and serum Cortisol 
and norepinephrine levels. Results: Serum Cortisol levels were significantly de- 
creased at 2 weeks compared to baseline (/j>=0.003). Plasma norepinephrine levels 
at 4 weeks were significantly decreased compared to baseline (/?=0.010). Heart 
rate, using the power spectra, increased significantly after 2 weeks compared to 
baseline. Of autonomic nerve conduction measures, latency was significantly in- 
creased at 2 and 4 weeks compared to baseline (p=0.023, 0.012), and amplitude 
was significantly decreased at 4 weeks compared to baseline (/?=0.008). There 
were no serious adverse events such as burns or other major complications. Con- 
clusion: The results of this study suggest that heat and massage applications pro- 
vide relaxation to the autonomic nervous system without serious adverse events. 
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Stress causes mental, emotional, and physical responses in humans. When stressed, 
an organism first experiences arousal of the sympathetic nervous system, followed 
by activation of the musculoskeletal, cardiovascular, and endocrine systems. The 
outcome of this process is a series of non-specific psycho-physiologic changes, 
such as an increase of certain neurotransmitters in the bloodstream, and results in 
adverse effects associated with the negative aspect of stress. Psychologically, high- 
ly differentiated synthetic behavior often occurs as a result of the activation of cog- 
nitive and emotional processes. Accordingly, when an individual is stressed, they 
may feel emotions such as fear and anxiety together with physiological responses. 
If this stimulus is repeated, anxiety may become learned or conditioned, and the 
response may become generalized. At present, the precise mechanisms of the in- 
teraction between stress and the responses of various systems have not been fully 
elucidated, nevertheless, it is generally recognized that systemic responses may in- 
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fluence the maladaptive behavior of an organism in re- 
sponse to prolonged stress. For this reason, examining ther- 
apies that modulate stress responses in organ systems such 
as the musculoskeletal, nervous, and endocrine systems 
may be helpful in psychiatric illnesses such as depression 
and generalized anxiety disorder. 1 

Heat application is thought to enhance metabolism and fa- 
cilitate circulation by dilating blood vessels and improving 
the functions of enzymes. These actions lead to increased 
catabolism, excretion of lactic acid, free fatty acids, and sub- 
cutaneous fat, and the removal of uric acid and other acidic 
waste products from muscle cells. By these mechanisms, 
heat therapy is thought to reduce fatigue and signs of aging, 
and to produce an analgesic effect. Heat is clinically effec- 
tive for the management of various injuries and to relieve 
back, shoulder, and other types of musculoskeletal pain. 2 

To decrease allodynia and stress in clinical practice, mas- 
sage is sometimes used as alternative management. When 
healthy volunteer or patients have been given massage, 
change in autonomic nervous system such as heart rate 
variability (HRV), blood pressure and respiratory rate have 
been observed. These results suggest that massage influenc- 
es the autonomic nervous system and changes stress reac- 
tions. Despite recent researches in this field, however, only 
a few models of the effect of massage on stress responses 
have been reported. 3 One hypothesis is that moderate mas- 
sage can increase response in the autonomic nervous sys- 
tem by activating sensory receptors, 4 and another is that 
touch massage influences to the release of oxytocin. 5 

The aim of this study was to evaluate the effects of heat 
and massage application on autonomic response in relation 
to stress, as measured by HRV, sympathetic skin response 
(SSR), and serum hormone (Cortisol and norepinephrine) 
levels in healthy adults. 



MATERIALS AND METHODS 



Study design 

This study has been approved by Institutional Review Board 
at the Yonsei University Wonju Christian Hospital. Ethical 
considerations were addressed by explaining the research to 
the participants, and obtaining their written consent. In re- 
sponse to advertisements on bulletin boards and local hos- 
pital newspapers and to verbal requests, 151 volunteers ap- 
plied for participation. At each session, if the participants 
refused to participate due to unwillingness to continue, heat 



and massage intervention was not carried out. Participants 
were excluded from the study if they have history of acute 
illness, malignancy, pregnancy, muscle disease, spinal pos- 
terior arch defect (e.g. spina bifida), previous operative his- 
tory of lumbar and thoracic back region, and cognitive dys- 
function. Demographic data were recorded for all subjects. 
In order to evaluate exclusion criteria, all subjects also un- 
derwent physical examination, electrocardiogram, radiogra- 
phy of the lumbar region, computed tomography of the 
lumbar region, and electromyogram (EMG). 

The physiologic effects of heat and massage application 
vary, depending on tissue temperature, the duration of ap- 
plication, the rate of tissue temperature increase, the size of 
the application area, and location. To minimize this effect, 
the study was carried out when temperature was maintained 
50-55°C for 40 minutes, and automatically applied on T9- 
L3 spinal ganglia, in order to provide similar pressure on 
the same locations for each subject and session. To main- 
tain constant temperature, temperature was constantly mon- 
itored with thermometer in the machine. To maintain con- 
stant massage area, the desired spinous process level was 
determined using the L5-S1 intervertebral space as an ana- 
tomical landmark. Furthermore, a line connecting the supe- 
rior iliac crests located the L4 vertebra. Heat and massage 
applications were performed with Ceragem M3500 (CER- 
AGEM Co. Ltd., Cheonan, Korea). The ergonomically de- 
signed guiding rail is attached to the horizontal mobile parts 
of the bed. This enables the internal projector to allow the 
massage effect. Heat was generated by the rolling jades 
(Fig. 1). It was carried out for 40 minutes daily, 5 days a 
week for 2 weeks, setting the temperature to 50-55°C at the 
hospital care room. The following variables were recorded: 
SSR, HRV, and serum Cortisol and plasma norepinephrine 
levels. Physical examination findings and vital signs were 
recorded at the beginning of the study and at 1 week after ini- 
tial application. After 2 weeks of application, subjects again 
underwent physical examination, and HRV, SSR, serum Cor- 
tisol, and plasma norepinephrine were measured. Two weeks 
after the cessation of application, and 4 weeks after the be- 
ginning of the study, these measures were repeated and re- 
corded (Fig. 2). Intertester reliability on all dependent mea- 
sures was established between clinical evaluators before the 
initiation of the study. 

Monitoring heart rate variability 

The response in the sympathetic and parasympathetic ner- 
vous system was estimated based on analysis of HRV, 
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Fig. 1. The picture of treatment device (Ceragem M3500). The internal pro- 
jector moves along the natural back line. 

HRV Cortisol, HRV Cortisol, HRV Cortisol, 

Norepinephrine Norepinephrine Norepinephrine 

SSR SSR SSR 
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Fig. 2. Schematic model, describing the study procedure. 

which refers to the beat-to-beat fluctuations in heart rate. 
The parasympathetic nervous system can quickly and finely 
adjust the time instant for the next heart beat, whereas sym- 
pathetic nervous system is a slower system for regulation of 
heart rate. This difference in activities creates variability 
within different frequency domains. The high frequency re- 
gion (0.15-0.40 Hz) is mainly mediated by parasympathetic 
activity, whereas the low frequency region (0.04-0.15 Hz), 
standard deviation of RR interval (SDNN) and total power 
(TP) are mainly mediated by sympathetic nervous system in 
relation to stress. HRV was analyzed by power spectrum 
analysis of interbeat intervals. RR intervals were converted 
to a time series by cubic spline interpolation, followed by 
resampling at 2.4 Hz. HRV was determined using an SA- 
2000E (MEDICORE Co. Ltd., Seongnam, Korea). The fol- 
lowing HRV indices were calculated: standard deviation of 
all normal RR intervals: the TP: the normalized low fre- 
quency power (LFP). 

Monitoring stress hormones (serum Cortisol and plasma 
norepinephrine) 

To obtain quantitative measurements of autonomic nervous 
system response, Cortisol and norepinephrine levels were 
measured from a single 10 mL blood sample, collected into 
serum and plasma blood tubes. Participants were instructed 
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to fast 2 hours prior to test, and blood sample were collected 
at same time of the day to minimize circadian and diurnal 
fluctuations in the biomarkers. Plasma was mixed and spun, 
and serum was allowed to clot for 30 minutes, and they were 
then centrifuged for 10 minutes. The quantitative measure- 
ment of Cortisol and norepinephrine was performed using a 
commercially available competitive binding enzymimmu- 
noassay according to the manufacturer's instructions. Data 
were collected three times at the beginning of study (pre- 
test), 2 weeks after application, and 4 weeks after applica- 
tion (2 weeks after cessation of application) (Fig. 2). All 
analyses were conducted in the Clinical Chemical Labora- 
tory at the University Hospital. 

Sympathetic skin response test 

SSR test has proven to be useful in examining the function 
of sympathetic activities. As a reflex response, an abnormal 
SSR might be due to the impairment of afferent and effer- 
ent peripheral nervous system, and psychological dysfunc- 
tion. Increased sympathetic activities such as high amplitude 
and decreased latency in pathologically stressful subjects 
have been observed. 

Using a Viking IV (CareFusion, San Diego, CA, USA) 
instrument, SSR was monitored as an additional indicator 
of the change of the sympathetic nervous system reflecting 
stress. In order to analyze the neurocutaneous response in 
the frequency range, the filter setting was adjusted (0.1-20 
Hz). To study the SSR, the subject was asked to relax and 
breath regularly for a short period of time, lying on a bed in 
a quiet and dimly lit room with an ambient temperature be- 
tween 22 and 24°C The skin underlying the recording sites 
was cleaned carefully with a dry cloth immediately prior to 
placing the surface electrodes, which were attached with 
sticking plaster. A round surface electrode 1 cm in diameter 
was attached 2.0 cm from the medial border of the left 
palm, and the reference electrode was attached at the oppo- 
site point on the dorsal side of the hand. Nerve stimulation 
was applied at maximum intensity on the proximal part of 
the volar wrist in the right median nerve distribution. 

Statistical analysis 

Repeated ANOVA test was used to analyze data between 
evaluation times. For data analysis before and after heat and 
massage application, paired-t tests were used. Analysis was 
performed using the SAS statistical software (Version 9. 1. 
3) and statistical significance was accepted for /?-values less 
than 0.05. 
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RESULTS 



General characteristics of subjects 

Among 151 initial subjects, three subjects were excluded 
because of arrhythmias found on the initial electrocardio- 
gram (one subject with first-degree artrioventricular block 
and two subjects with premature ventricular complexes) 
and nine subjects dropped out during the study, yielding a 
total of 139 subjects who completed the study. Image anal- 
ysis and EMG of the lumbar region were performed on all 
subjects for anatomical electrophysiological evaluation that 
include exclusion criteria, e.g. spina bifida, spinal tumor and 
muscular disease. There were no participants that had exclu- 
sion criteria; however, 54.3% had the following findings: in- 
Table 1. General Caracteristics of the Subjects (n=1 39) 



Variables 


Number of subjects (%) 


Sex 


M/F 


8(5.8)/131 (94.2) 


Age (yrs) 


21-30 


1 (0.7) 


31-40 


27 (19.4) 


41-50 


55 (39.5) 


51-60 


33 (23.8) 


61-70 


15 (10.8) 


71-80 


8(5.8) 


CT finding 


Normal 


55 (39.5) 


Disc protrusion 


43 (30.9) 


Spinal stenosis 


7(5.1) 


DS 


23 (16.6) 


Others 


11(7.9) 


EMG finding 


LR 


14(10.1) 


Normal 


125 (89.9) 



F, female; M, male; DS, degenerative spondylosis; EMG, electromyogram; 
LR, lumbar radiculopathy. 
Values are numbers (percentages). 



tervertebral disc protrusion in 43 subjects (30.9%), degenera- 
tive spondylosis in 23 subjects (16.6%), spinal stenosis in 
seven subjects (5.1%), spondylolysis in seven subjects 
(5.1%), intervertebral disc calcification in four subjects 
(2.9%), and lumbar radiculopathy in fourteen subjects 
(10.1%) (Table 1). 

Analysis of variability 

Analysis of serum Cortisol 

The average serum Cortisol levels of the subjects were as 
follows: 9.54±2.01 pg/mL at baseline, 6.92±1.62 pg/mL af- 
ter 2 weeks of therapy, and 8.29±2.11 pg/mL at 4 weeks af- 
ter initiating the study. Serum Cortisol levels after 2 weeks 
of treatment were significantly decreased compared with 
the baseline (/?=0.003). Serum Cortisol measured 4 weeks 
after beginning the study tended to be decreased compared 
with the baseline, but these results were not statistically sig- 
nificant (Fig. 3). 

Analysis of plasma norepinephrine 
Average plasma norepinephrine levels of the subjects were 
as follows: 190.5±57.9 ug/dL at baseline, 170.9±93.3 ug/dL 
after 2 weeks of therapy, and 132.8±41.1 (J.g/dL at 4 weeks 
after beginning the study. These levels were significantly 
decreased at the 4-week time point compared with the base- 
line (p=0.010) (Fig. 3). 

Analysis of power spectrum of heart rate change 
In the analysis of the power spectrum of the heart rate change 
index after 2 weeks of heat and massage application, TP and 
SDNN were significantly increased (p=0.014, 0.021), but 
there was no significant difference in the degree of change, 
SDNN, and TP after 4 weeks compared with the baseline. 
Normalized LFP tended to generally decrease, but the find- 
ings were not statistically significant (Fig. 4). 
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Sympathetic skin response 

In the analysis of SSR, the latencies at 2 and 4 weeks post- 
treatment were 1,405.4±150.4 msec and 1,463.4±138.9 msec, 
respectively. Compared with the baseline of 1,231.0±221.3 
msec, the latencies after 2 and 4 weeks were significantly 
increased (p=0.023, 0.012). Amplitude was 1,873.7±780.9 
uV after 2 weeks and 1,598.2±936.2 uV after 4 weeks post- 
treatment, and was significantly decreased compared with 
the baseline of 2,108.4±1193.1 uV (p=0.008) (Fig. 5). 

Adverse events 

Adverse events were evaluated at 2 and 4 weeks, and there 
were no serious adverse events such as bum or other major 
complications. Ten subjects (7.2%) complained of a tempo- 
rary increase in low back pain. 



DISCUSSION 



The aim of this study was to evaluate the effects of heat and 
massage application on autonomic responses, especially 
sympathetic activity, in relation to stress. The main findings 
in this study were that, after 2 weeks of heat and massage 
application, there was significant decrease in serum Cortisol 
level, and there was significant decrease in serum norepi- 
nephrine level after 4 weeks, indicating an overall down reg- 
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ulation of sympathetic activity. In measurements of HRV us- 
ing the power spectra, there was significant increase in TP 
and standard deviation, mainly mediated by sympathetic 
nervous system. Of autonomic nerve conduction measures, 
latency was significantly increased at 2 and 4 weeks, and 
amplitude was significantly decreased at 4 weeks. These 
findings would normally be interpreted as decreased down 
regulation of sympathetic activity, indicating reduced stress 
response. Ten subjects (7.2%) complained of a temporary 
increase in low back pain. The cause to induce temporary 
increased low back pain among our study population is 
likely multifactorial. First, although temperature and pres- 
sure were maintained to induce heat and massage effect, it 
is possible to induce tissue damage and consequent noci- 
ceptor activation. Second, it has been suggested that pre-ex- 
isting low back discomfort could be extensive. 

Activation of the autonomic nervous system, including 
hemodynamic, electrodermal, and hypothalamic-pituitary- 
adrenal system, influences metabolism. The sympathetic 
nervous system and stress contribute to catabolic activity and 
the consumption of energy, whereas the parasympathetic 
nervous system works for anabolic activity and the storage 
of energy. Stress can, therefore, cause catabolic activity. Our 
findings support the hypothesis that HRV, serum Cortisol 
and norepinephrine levels, and SSR would be lowered after 
heat and massage application. 
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The autonomic nervous system is responsible for main- 
taining homeostasis against external conditions by control- 
ling the activities of the viscera, blood vessels, and secretory 
glands. 6 As the incidences of certain diseases such as diabe- 
tes, affecting the vascular system, are increasing, the aware- 
ness of their effects on the autonomic nervous system has 
been heightened. Various evaluation methods have been de- 
veloped to examine the function of the autonomic nervous 
system, including measuring cardiovascular parameters re- 
lated to posture and isometric exercises, measuring body sur- 
face temperature with infrared ray photography, evaluating 
the mobility of the gastrointestinal tract using certain iso- 
topes, analyzing blood levels of autonomic neurotransmit- 
ters, testing electrophysiological parameters, examining the 
autonomic nervous system response to medication, quantita- 
tively measuring sweating, and evaluating pupils and puden- 
dal nerve responses. 7 9 However, most of these evaluation 
methods have limited clinical applications because of diffi- 
culty in quantifying the results, lack of ability to reproduce 
the results, or they are very invasive. 10 Heart rate is constant- 
ly changing to maintain homeostasis, and is determined by 
autonomic nervous system stimulation of the S-Anode and 
spontaneous excitation of the S-A node. 1112 Since heart rate 
is controlled by the antagonistic work of sympathetic and 
parasympathetic nervous systems, the analysis of heart rate 
change can reflect the balance of each component of the au- 
tonomic nervous system. 61318 This analysis yields three 
points on the power spectrum: an ultra low frequency com- 
ponent of less than 0.05 Hz related to thermoregulation and 
the renin-angiotensin system, a low frequency component 
at around 0. 1 Hz, and a high frequency component at around 
0.25 Hz. 10 - 12 - 17 - 19 ' 24 

In this study, analysis of the SDNN and TP indices were 
significantly increased after 2 weeks of heat and massage 
application, compared with the baseline, however, they 
tended to decrease 2 weeks after heat and massage applica- 
tion stopped, that is, 4 weeks after the study began. LFP was 
found to be generally decreased, but the findings were not 
statistically significant, implying that heat and massage ap- 
plications lead to up-regulation of both autonomic activity. 
This notion is supported by the view that sympathetic and 
parasympathetic branches do not always act reciprocally, 
but may also act synergistically and complementarily. 

The mechanism of the SSR has not yet been fully de- 
fined, but it has been demonstrated that electrical stimula- 
tion of peripheral nerves is transmitted to the lateral col- 
umns of the spinal cord through group II, and III fibers, that 



sweating is caused by transmission of impulses to the neu- 
ro-sweat gland junction through pre-ganglion and post-gan- 
glion small diameter centrifugal fibers, and that temporary 
skin electric potential changes due to sweating can be re- 
corded and are known as the SSR. 13 - 2526 Since Knezevic 
and Bajade 27 outlined a method to record this skin response 
non-invasively, many researchers have been studying it. 
Aisen and Stallman 28 reported that spinal cord injuries in- 
fluence the SSR after examining it in patients with paraple- 
gia, and Brown and Wang 29 discovered that sensory motor 
areas, hypothalamus, and brainstem reticular system all ex- 
ert an excitative effect on the SSR. The connections within 
the central nervous system for the provocation of the SSR 
consist of a multisynapse pathway, but the connection be- 
tween the hypothalamus and the motor neurons of the spi- 
nal cord is unclear. The descending autonomic pathway in 
the brain stem and spinal cord is not limited to one area but 
is widely spread throughout. 30 - 31 In 1993, Cheong, et al. re- 
ported that the latency of the SSR in spinal cord injury pa- 
tients was longer and its amplitude was greater than that of 
normal controls. There was also a report that, in the case of 
patients with cerebrovascular accident, if the nerve stimulus 
area and the record area were different between right and 
left, the latency was longer than when they were same. 32 
Furthermore, in patients with anxiety and multiple sclero- 
sis, 94.2% of SSR were found to be abnormal, and the ratio 
of abnormal findings was higher than that of any other 
evoked potential test. 30 Exposure to continuous stress causes 
functional changes in the autonomic nervous system, and 
the SSR intensity in electric nerve conduction tests is gener- 
ally recognized to increase such as increased amplitude and 
decreased latency. 33 In SSR measurements in this study, the 
latency was found to be significantly increased after 2 and 4 
weeks of heat and massage application, and the amplitude 
was significantly decreased compared with the baseline af- 
ter 4 weeks of application, indicating overall decrease of 
sympathetic activity. 

Catecholamine release in the bloodstream is increased 
under a variety of conditions, including stresses, hypogly- 
cemia, hypovolemia secondary to bleeding, and hypoxemia. 
Furthermore, it is very sensitive to changes in the blood 
sugar level, posture, and activity level, and also shows epi- 
sodic daily changes. Cortisol is secreted from the adrenal 
cortex, in response to ACTH release by the pituitary gland, 
and more than 90% of it in the blood are conjugated to pro- 
teins and excreted in urine after being metabolized in the 
adrenal gland, liver, and kidneys. Together with aldoste- 
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rone, Cortisol is an important hormone in the adrenal cortex 
and an important part of the stress response. 31 Norepineph- 
rine secreted by sympathetic nerve endings binds to recep- 
tors on neighboring organs and influences their function, 
and then undergoes re-uptake at the sympathetic nerve end- 
ings. Epinephrine and norepinephrine secreted by the adre- 
nal medulla enter the bloodstream and exert their effects on 
organs farther away. These hormones work by stimulating 
the a and P receptors, each having a specific effect at each 
receptor type. Adrenal medulla secretes epinephrine and 
norepinephrine at the ratio of 9 : 1, measured in the renal 
vein. In the capillaries, this ratio decreases to about 1 : 4. 
This is due to the fact that epinephrine is made only in the 
adrenal medulla and the clearance of epinephrine and nor- 
epinephrine is almost the same. Epinephrine secretion can, 
therefore, be used as a reflex of sympathetic nerve function, 
as well as adrenal medulla function. In clinical diagnostic 
testing, blood catecholamine levels can be used to diagnose 
pheochromocytoma and neuroblastoma. In addition, cate- 
cholamine increases can be seen in essential hypertension, 
renal hypertension, malignant hypertension, and congestive 
heart failure, but catecholamine monitoring in these diseas- 
es is not diagnostically useful. Under stress, catecholamine 
levels can also be increased. 34 In this study, plasma norepi- 
nephrine after 2 weeks of heat and massage applications 
was found to decrease at 2 and 4 weeks compared with the 
baseline. Serum Cortisol levels at 2 weeks and plasma nor- 
epinephrine levels at 4 weeks were also significantly de- 
creased. These findings indicate decreased activity of the 
Hypothalamic-Pituitary-Adrenal axis. These implicate that 
heat and massage application contribute to decrease cata- 
bolic activity. 

In conclusion, we objectively evaluated that the effect of 
heat and massage application provide relaxation in the au- 
tonomic nervous system through comparison analysis of 
data obtained before and after heat and massage applica- 
tion. Despite the impressive results of heat and massage ap- 
plication found in this study, cautious interpretation of 
study findings is warranted because of some limitations, 
such as experimental design without a control condition, 
study design without a self-report measure of stress to have 
clinical application. The majority of the sample was female 
so there was some limitation to generalize to male popula- 
tion. These results reflect the effects of heat with mechani- 
cal massage but cannot be compared directly to results of 
manually-derived massage. Given limitations noted in this 
study, future studies with blinding and randomized alloca- 



tion of study subjects, a controlled group, and the effect of 
heat vs. massage isolated study are important to validate the 
effectiveness of heat and massage application. 
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